Pairing symmetry and vortex zero-mode for superconducting Dirac fermions 
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We study vortex zero-energy bound states in presence of pairing between low-energy Dirac 
fermions on the surface of a topological insulator. The pairing symmetries considered include the 
s-wave, p-wave, and, in particular, the mixed-parity symmetry, which arises in absence of the inver- 
sion symmetry on the surface. The zero-mode is analyzed within the generalized Jackiw-Rossi-Dirac 
Hamiltonian that contains a momentum-dependent mass-term, and includes the effects of the elec- 
tromagnetic gauge field and the Zeeman coupling as well. At a finite chemical potential, as long 
as the spectrum without the vortex is fully gapped, the presence of a single Fermi surface with a 
definite helicity always leads to one Majorana zero-mode, in which both electron's spin projections 
participate. In particular, the critical effects of the Zeeman coupling on the zero-mode are discussed. 

PACS numbers: 03.75. Lm,74.20.Rp,74.45.+c 



I. INTRODUCTION 

The zero-energy Majorana modes have recently be- 
come a topic of many investigations in condensed mat- 
ter physics, perhaps due to the concomitant non-Abehan 
statistics [l| and potential applications to fault-tolerant 
quantum computation[2|. In a weakly-coupled triplet 
p-wave superconductor, for example, one such zero- 
mode was found in the core of a half-vortex, and two 
pairs of such vortices are required to perform the non- 
Abelian braiding operations [3]. The vortex bound states 
were, of course, already studied a long ago by Car- 
oli, de Gennes, and Matricon in a s-wave weak-coupling 
superconductor,!^] but a pair of true zero-modes was dis- 
covered only much later, by Kopnin and Salomaa in a 
triplet superfluid 5] . Both studies were performed within 
the framework of Bogoliubov-de Gennes (BdG) equations 
which describe the low-energy electron-hole excitations 
around the degenerate Fermi surfaces in the presence of 
superconducting instability. The degenerate Fermi sur- 
faces can be split, and the superconducting order pa- 
rameter can contain both the spin singlet and triplet 
components when the inversion symmetry is broken[6J. 
For example, the corresponding vortex bound states in a 
noncentrosymmetric superconductor, ;Jl] or the Andreev 
bound states on the domain wallQ have also been stud- 
ied. The existence of such zero-modes has been shown 
to depend on a Z2 topological invariant associated the 
mixed order parameter p3l4ll| ■ 

The metallic surface of a topological insulator (TI) 
resembles the noncentrosymmetric material, as in both 
cases the Fermi surfaces are chiral, i.e. the projection 
of the spin is locked to the direction of the momen- 
tum. However, an important difference is that the surface 
of the topological insulator has a only one such Fermi 
surface [iJI- Recently, an interesting system consist- 
ing of this peculiar surface subjected to the supercon- 
ducting proximity effect has been proposed to bind one 
Majorana zero-mode inside a s-wave vortex core of unit 
vorticityjlSJ- This zero- mode is nothing but the mid-gap 
state of the Jackiw-Rossi-Dirac Hamiltonian jl4| (JRD), 



studied before in the context of charge fractionalization 
in polyacetylene [15| and graphene.flGl [l7| with its exis- 
tence assured by an index theorem[l8|. In addition, the 
zero-modes of the JRD Hamiltonian have also been stud- 
ied recently in the context of the insulating [la, [I^, [20| 
and superconducting states on the graphene's honeycomb 
latticei2l|-|25J. 

In this paper we study the existence and the form of 
the zero-modes of a class of the general JRD Hamilto- 
nians describing paired Dirac fermions in presence of a 
superconducting vortex in the spin-triplet p-wave state, 
and in the mixed singlet-triplet state, when the inver- 
sion symmetry is absent. We pay special attention to 
the effects of a finite chemical potential[2j, [2J| and the 
magnetic fields on the zero mode, which are all present in 
real systems. The problem under study is also relevant 
to graphene, where different symmetries of the super- 
conducting phases have also been considered |26l - [28| . In 
graphene in particular, other exotic manifestations of the 
physics of zero-modes become possible due to their rich 
internal structure, consisting of spin, valley, sublattice, 
and the Nambu degrees of freedom [1^, [l^, [2^ [29|. It is 
therefore of importance to understand the general condi- 
tions for the appearance of the zero-energy states. Quite 
generally, we find that the zero-modes in the supercon- 
ducting vortex may exist only for weak enough Zeeman 
coupling to the magnetic field. This is in accord with our 
previous study limited to the the vortex in the s-wave su- 
perconducting order parameter [30l|. 

This paper has the following organization. We first in- 
vestigate the spectrum associated with the spatially ho- 
mogeneous Hamiltonian without vortex for the singlet 
s-wave, triplet p-wave and mixed-symmetry cases. The 
conditions for opening of a gap in all three cases are dis- 
cussed. We then proceed to the equations for the zero- 
mode associated with different pairing symmetries, with 
the Zeeman term and the chemical potential included. In 
the sec. IV, the effects of the Zeeman field and the rela- 
tion of the zero-mode with the antilinear operator that 
provides the reflection symmetry of the spectrum [30| in 
all symmetry cases are discussed. Finally, a brief sum- 




FIG. 1: A single two-dimensional Fermi surface associated 
with the Dirac Hamiltonian, Ho{k) = VF{i ■ k x a) — fi, on 
the metallic surface of a TI. z, the out-of-plane direction, is 
the spin quantization axis. The dark (dashes) arrows indicate 
the directions of the spin associated with the momentum state 
k on the Fermi surface of positive (negative) helicity when fj, 
is positive (negative). 



mary is given. 



II. TRANSLATION-INVARIANT MASSIVE 
DIRAC HAMILTONIAN 



For later convenience, we note that H{k) can also be writ- 
ten in terms of the cross-product representation, namely 
H(k) = {vpnik) ■ a - ^i)T3 + M{k)T+ -f M^k)T^, where 
CTj's act on the spin and r^'s on the Nambu indices. The 
vector n(fc) = (— fcy, kx) hes on the x-y plane, as in Fig. [T] 
The chemical potential jj, enters along with the number 
operator N, which can be identified as ra in this repre- 
sentation. As we discuss shortly, while they of course 
never commute with the number operator t^, the mass 
terms may commute with the diagonal, kinetic energy 
terms. These cases will lead to a peculiar spectrum that 
is not possible in standard superconductors which have 
the quadratic energy-momentum dispersion in its normal 
state. 

First, let us discuss different pairing symmetries for the 
mass term, in the translation-invariant situation without 
the vortex, and the resulting spectrum E{k). 

The simplest case corresponds to the singlet s-wave 
pairing in which A is constant and ua=0. The spectrum 
may be obtained most easily by squaring the Hamiltonian 
and using the anticommuting properties of the Pauli ma- 
trices; this gives [H{k)Y = {vpn{k) ■ a — ^)^ -|- M^, and 
therefore E{}i) = ±^{vFk±nY + A2.[ll, ^ Next we 

consider the triplet pairing, with the vector d coinciding 
with n, in which the mass term M = vau ■ ff still anti- 
commutes with the kinetic energy. The same trick can 
be applied again, to obtain the spectrum. 



In this section, we consider the Dirac electrons on the 
metallic surface of a TI and their pairing due to the prox- 
imity to a superconducting material. The metallic sur- 
face of a TI has a single Fermi surface as shown in Fig. [TJ 
and the spin-dependent spectrum can be described by a 
simple 2x2 Hamiltonian, Ho{k) = vp^ ■ k x a — fi^ where 
fj, is the chemical potential and vp is the Fermi velocity. 
The quantization axis for spin is parallel with the normal 
of the surface z. As usual, the identity matrix and the 
three Pauli matrices, {/2,(?}, define a basis in the two- 
dimensional matrix space. One should note that Ho{k) is 
invariant under the point group C^y in which the inver- 
sion symmetry is absent [31]. At the mean- field level, the 
general pairing potential Vp that respects the absence of 
inversion symmetry and takes different coupling channels 
(spin singlet and triplet) into account needs to be be ex- 
pressed as a matrix. [Q] To this purpose we introduce the 
general mass matrix M{k) — AI2 + VAd{k) ■ a, so that 
the pairing potential becomes Vp{k) = iM{k)(jy. Here A 
stands for the singlet pairing gap while the velocity ua 
is associated with the triplet pairing. 

In the 4-dimensional spin-Nambu space, the BdG 
Hamiltonian 



H = Y.k%H{k)^., with 



fc^' 



J. *} , then assumes the form. 



H{k) 



Hoik) 

Mt(fc) 



M{k) 
-Hoik) 



(1) 



Eik) = ±^ivFkTti)^ + ivAk)\ 



(2) 



which is gapped as long as /i is finite. At /i = 0, the 
spectrum is gapless, and in fact reduces to the Dirac 
cone, with an increased velocity. A somewhat unusual 
situation arises when d = ikx,ky), for which the mass 
term commutes with the kinetic energy: when n_L(i, 
[MT±,n ■ i7T3]=0. It can be shown that the spectrum is 
then given by Eik) = -iivpk± \J [j? -j- (wAfc)^- In this case 
the spectrum is still gapless, and it is the chemical poten- 
tial that becomes modified and momentum-dependent, 
loosely speaking [33|. 

In the absence of the inversion symmetry the general 
"mixed" mass term, with both A and wa finite, is also 
allowed. Let us assume the triplet component to be given 
by d = n. By squaring the Hamiltonian, one can obtain, 
\H(ky^ = ivFfi-(T—^)^ + iAs+VAn-(T)^, and consequently 



Eik) - ±v/(«FfcTM)' + (A±i'Afc)2, 



(3) 



where a gapped spectrum exists only for the positive he- 
licity branch when A — v/^kp = A+/2, with the Fermi 
wavevector kp = [ijvp. 



III. EQUATIONS FOR ZERO-ENERGY 
QUASI-PARTICLE STATE 

In this section we shall consider the inhomogeneity in 
the mass term M generated by an isolated vortex in 
which the magnetic field applied along z is threading 
through its core. The superconducting phase changes by 
27tQ on circling once around the core. Here Q labels the 
winding number of the vortex. In the context of super- 
conductivity, the mass term here is in general a function 
of both the center-of-mass position r and relative momen- 
tum k of the Cooper pair, namely, M = M{r,k). More 
explicitly, in the polar coordinate r — (r cos (j), r sin (f) , 
the vortex enters the mass via Af (r, k) = e^'^'^M{k). In 
the following subsections, we shall employ the BdG equa- 
tion to investigate the zero-energy bound state in the 
vortex core for various pairing symmetries. In particu- 
lar, we also consider the effects of an external Zeeman 
field Hi = hY^-j:'^*-.[a-iTo)^ j:, in which tq stands for the 
identity matrix acting on the Nambu space. 



A. Singlet s-wave pairing 

The mass generated by the vortex of singlet s-wave 
pairing potential with the winding number of unity can 
be written as M = e*'^|A(r)|, in which the magnitude 
increases from zero at the core and reaches a positive 
value of Aoo when the distance r is sufficiently large 
than the superconducting coherence length. In terms 
of the Dirac fermion ^ we defined, the quasi-particle 
states which diagonalize the BdG Hamiltonian and give 
exactly zero-energy have the corresponding coefficients 
u = {ui{r),U2{r)} for the electron sector and the cor- 
responding components v for the hole sector satisfy the 
following equations. 



-li + vpi- {k-eA) xa + ha^ u + e^-^IAWIv = 0,(4) 
e-''^|A(r)|u+ L- Wi^z- {k + el) x a + ha^ v = ,(5) 



where ?i=:c=l is taken for simplicity. The momentum 
operators, ki = —idi, are coupled to the gauge field A — 
(— y, x)B{r)/2, where B varies on a length scale of pene- 
tration length Ip inside a superconductor. Using the iden- 
tity {dx ± idy)—e'^'"^{dr ± {i/r)d^), and the angular de- 
composition, [u, v]^=[Mi(r),e*'^U2(r),e~*'^Ui(r),i)2(r)]'^, 
one can decouple the set of four first-order differential 
equations into two blocks, namely. 



VF{dr + ---B) 

r 2 



lAI 



/„ + (m - h)fb = , (6) 



VFidr + -B)-\A\ /6-(/i + /l)/a=0, (7) 



r 2 



ga + ifi- h)gb = , (8) 



VF{dr + 7:B) + \A\ .gb-(M + %a = 0, (9) 



where the definitions U2 ± vi=faiga) and mi =F V2—fb{gb) 
are made. The explicit dependence of magnetic field B 
and pairing potential A on the radial coordinate r is not 
shown to simplify the notation. In fact, it will be seen in 
the following asymptotic analysis that the effects of cou- 
pling to the gauge field on the zero-mode are negligible. 
Far from the core the gap function A(r) — > Aoo and 
moreover, only the derivative terms are relevant in the 
kinetic energy. The terms related to B{r) can be dropped 
since B{r) decay faster than 1/r in the case of a isolated 
vortex 34] . Therefore, by defining the new variable x = 

h'^\r/vF, Eqs. 



following: 



([U and (O are reduced to the 



dx"^ 



W, 



Ux) 



= 0. 



(10) 



where the plus (negative) signs are indicative of that |/z| 
is larger (smaller) than |/i|, respectively. In both cases, as 
long as h? < fj?' + A^ , the two solutions for fa are both 
exploding at infinity, which, from Eq.(|6l), is also true for 
fb- Thus one concludes that ui — V2 and U2 — —vi, to 
ensure the normalizability. The remaining components 
can be determined then from Eq. ([S]) and © , which can 
be reduced to, 



(^*" 



, VU'^-h^ 



gaix) 



0, 



(11) 



where two decaying solutions are available under the 
same condition h? < ^"^ + A^ . 

Near the core where the argument x <^ 1, it is le- 
gitimate to keep only the singular and the derivative 
terms. This is because A(r) is vanishingly small there, 
and the inherent logarithmic divergence from _B(r) 34] is 
suppressed in the vector potential. Eqs. ([8]) and ([9]) can 
also be written as, 



cp Id 
dx^ X dx 



±l)gaix)=0, 



(12) 



where the plus sign still means that |/i| is larger than 
\h\. Now one can see that U2{x)=Ji{x), the first-order 
Bessel function, for small argument. Note that the other 
solution, TVi (x) , is forbidden because of its singularity of 
1/x for small argument. For |/i| < \h\, U2{x)=Ii{x), the 
modified Bessel function which is regular at small x. The 
other component ui{x) is not independent and is given 
by. 



ui{x) 



V% 



1, 



,^i^ + -)u2ix) , (13) 

II — h dx X 



for small argument x. If we think of the pair poten- 
tial A(r) as a step function which vanishes within some 



certain range and is a constant outside, the obtained 
solution from the second-order differential equation for 
cc <C 1 can be smoothly joined with the solution for x 3> 1 
at the boundary by matching the values of the function 
and its first derivative. This can be accomplished when 
h^ < IJp' + A^ since, given one function on the left-hand 
side of the boundary, there are two available decaying 
solutions on the right-hand side. For an arbitrary A(r) 
the solutions can be obtained by adiabatic deformation 
of the step-like pair potential that leaves the asymptotic 
behaviors intact. 

We conclude that, as long as h? < ^^-|-A^, there exists 
one zero-energy solution for the s-wave pairing vortex 
of winding number one, with the following radial wave 
function. 



/ Mr) 
U2{r) 
-W2(r) 
V ui(r) 



(14) 



where the electron components Ui and U2 for small argu- 
ment are given by, 



ui{x) 

U2{x) 



when h^ < ^^ . As for fj? 
Ji{x) are simply replaced with the corresponding modi- 
fied Bessel functions /o(a;) and /i(a;), respectively. For a 
large argument. 



-Mx) J 


(15) 


<h'<^^' + Al,, 


Jo{x) and 



ui{x) 

U2{x) 



VIa'^J' ' (^Acosx - Bsinx) 
— {A sin X + B cos x) 



(16) 

which is vahd only for h"^ < ^^ . For ^"^ < h"^ < fM^ + A^ , 
the trigonometric functions are replaced with the corre- 
sponding hyperbolic ones. The coefficients A and B are 
to match U2{x) at the boundary and an additional overall 
coefficient for normalization is not shown explicitly. 

A special case of above is when /i = /i = with the 
solution given by. 



be employed, and, consequently, ui ~ Ji and U2 ~ Jo 
near the core. Moreover, in the special case with vanish- 
ing /i and h, the Majorana zero-mode have similar form 
as in Eq. (|17p but with down-spin components only. 



B. Triplet p-wave pairing 



Now we consider the triplet order parameter specified 
by d = {—ky,kx), which is parallel to the spin direction 
n in the kinetic energy. With this Rashba type p-wave 
pairing, an isotropic superconducting gap is generated 
for nonzero fi according to Eq. © in contrast to the case 
with d J- n. The quasi-particle states with zero energy in 
the presence of vortex in the triplet order parameter with 
winding number of unity have to satisfy the following 
equations. 



WF(9r. -t- -)fa + {fi- h- VAdr)fb = , 

r 



Vpdrfb 



fj,+ h — VA{dr + -) 



/a = 0. 



1, 



vpidr + -)ga + ifj--h + VAdr)gb = , 
r 



VpdrC, 



^1 + h + VA{dr + -) 

r 



9a = 



(18) 
(19) 
(20) 
(21) 



where identical angular dependence has been employed, 
and fa{b) and ffo(b) have the same definitions as before. 
VA{r) is also a radial function. When ^ ^ h = 0, it 
can readily be deduced that both fa and ga must be 
identically zero; otherwise the operator (dr + 1/r) will 
lead to the divergence of ^ 1/r in the solution near the 
origin. Consequently, the zero-energy state must satisfy 
drfb — drgb = for all r, which leads to the conclusion 
that the zero-energy bound state does not exist when 
/i = 0. For general fi and h, on the other hand, we may 
proceed in the same manner as in previous subsection. 
Near the vortex core the analysis is identical since the 
gap velocity va is vanishingly small there. Far from the 
core, it is useful to rewrite Eq. (IT5|) and P^ in terms of 
fa ± ifb, as. 



ui{r) 
U2{r) 







^-^/'■|A(r')|dr' 



(17) 



which is exact for all range of r. Comparing with 
Eqs. ([T5|) and p^ . one should note that the ambi- 
guity arising at the zero chemical potential is removed 
by the fact that Ji is vanishing for zero argument. Be- 
sides, the zero-mode in presence of a finite chemical po- 
tential appears with both spin projections in contrast 
to the zero-mode in Eq. p7)) where only the up-spin 
components show up[l3|. For the case with an anti- 
vortex, M ~ |A(r)|e~*'^, different angular decomposi- 
tion, [u, v]^=[e~''^Mi(r),M2(f),wi(r),e**U2(r)]'^, has to 



[{VF 


+ iVA)dr 


- «M] {fa + 


^fb)' 


-lh{f^ 


*/fc) = 


= 0, 


(22) 


[{VF 


- iVA)dr 


+ T-lA {fa - 


^fb) 


+ ^/^(/ + ifb) = 


= 


(23) 



which also can be recast as a second-order differen- 
tial equation. For fj?' — h? > 0, the desired equa- 
tion is {d'^/dy^ — 2\d/dy + l)(/a + ifb) = 0, where the 
variable y = -^/j/i^ — h'^\/(vp + v'^)r and the constant 
A = IJ'Va/ \/{vp + f^)|/i^ — h'^\. The solution then has a 
plane wave form, e*'^^, with k = —i\ ± \/l — X^. With- 
out loss of generality va may be assumed to be posi- 
tive. We then first consider /x to be positive: there are 
then two exploding solutions for fa + ifb at infinity when 
< A < 1, or equivalently /i-^ > (1 + v\/vp)h'^. For the 



case of A > 1, or /i^ is more close to /x^, the two solu- 
tions are still exploding but with k = —i{X =F VA^ — 1). 
On the other hand, for /i^ < h^, the differential equation 
becomes {(P /dy"^ — 2Xd/dy — l)(/a + ifb) = 0, with the 
asymptotic solution of wave vector k = — i(A±-\/A^ + 1). 
In this case one solution is always decaying, while the 
other is exploding. Thus we conclude that fa = /& = 
as long as fi^ > h^, and we take ui = V2 and U2 ~ —v\. 
The remaining task then is to obtain the Mi's from the 
equations for ga and g;,, which are identical to Eqs. (P^ 
and (I23|) . except that ua -^ —va and hence two decaying 
solutions at infinity are available as long as /i^ > h^. 

For /i negative, the above conclusion has to be reversed, 
so that ga = gt = when fp > h^, and therefore ui — 
—V2 and M2 — vi. The rest of analysis is the same as with 
positive /x. The complete wave function for zero-energy 
state has the radial part summarized as following: 



■ui(r) 

U2ir) 

-sgn{^i)u2{r) 

sgn(pi)Mi(r) 

and the asymptotic behaviors are given by 



(24) 



uiiy) 
My) 



{A-B)smVr^ 
{A + B) cos ^/^ 



A2y 



A2y 



-Wv 



(25) 



which has an oscillation on top of a decaying form when 
< A < 1. It will be monotonically decaying at infinity 
when A > 1, and the trigonometric functions need to 
be replaced with the corresponding hypergeometric ones. 
The coefficients are determined by matching the function 
M2 + iui to the solution near the core. 



C. Mixed s-wave and p-wave pairing 

Finally, we consider the combination of the problems in 
the previous two subsections. In absence of the inversion 
symmetry the kinetic energy in the normal state can have 
its chiral form with the pairing potential as the mixture 
in two order parameters of different parities. We will 
thus assume it to be a mixture of the singlet and triplet 
components, with the latter specified hy d = {—ky,kx)- 
For simplicity, we also neglect the Zeeman effect. Far 
from the vortex core, the corresponding equations for /'s 
are, similarly as before. 



{vpdr - Aoo)/a + (al^l 

(vpOr - Aoo)/b - {a\^i\ 



VAdr)fb = , 
VAdr) fa = , 



(26) 
(27) 



in which a specifies the sign of chemical potential, and 
we take a — 1 first. The above equations are identical to 
those for ga and gt, except that Aoo — > — Aoo and va — > 
—VA- It follows that at infinity both fa ± ifh are in the 



form of e^'*='=*«)'' with k = {a\n\vA + AooW_f)/(va + ^f) 
positive, and q — {a\^j,\vF — AooWA)/('yA + Vp)- It is 
useful to define the parameters 



A± == Aoo ± VAkp 



(28) 



where vpkp = \fj,\ is the Fermi momentum for nonzero 
chemical potential and the renormalized Fermi velocity 
vp = vf{1 + v'^/vp). The exploding solutions can then 
be written as ~ g^'^+r/vp^ which is similar to the case 
of s-wave pairing. Of course, the exploding solutions 
are forbidden, and therefore ui = V2 and U2 = —vi. It 
can be seen that the set of equations for ga and gt will 
lead to two asymptotic decaying solutions in the form 
g-A+r/tiFgiigr^ -^^liich in turn determines the zero- mode 
wavefunctions to be. 



ui{r) 
U2{r) 






Be-'i"- 



- /S.-]-r /vf 



(29) 

with the coefficients determined by the similar boundary 
conditions for matching U2 to its piece near the core. 

For a negative chemical potential, by similar reasoning, 
the set of equations associated with /'s gives two asymp- 
totic solutions as e^-''/^^e***'', which will be forbidden if 
A_ > 0. Then, the zero-mode is determined from the set 
of equations associated with the g's, which give similar 
solution as in Eq. (|29p . except that A+ — > A_. However, 
the zero-mode will have ui = —V2 and U2 = Vi in contrast 
to the previous case when A_ < 0, in which the solutions 
from the /'s were decaying and allowed, while those from 
g^s were forbidden. The wavefunction for /z < therefore 
has a general form, 



ui{r) 

U2{r) 

-sgn(A_)M2(r) 

sgn(A_)Mi(r) 



(30) 



where one can see that a vanishing A_ will remove the 
zero-mode. The corresponding ui and U2 are determined 
in a similar manner. 



IV. DISCUSSIONS AND CONCLUSIONS 

In the preceding sections, the existence of zero-energy 
bound states in the presence of various symmetries of 
the mass-gap is investigated in the context of supercon- 
ducting Dirac fermions on the surface of a TI. It may 
also be useful to write the BdG Hamiltonian in pres- 
ence of the vortex in terms of the 4-dimensional Dirac 
F-matrices. The s-wave vortex, for example, is described 
by the Hamiltonian H = Fipi -l-r2P2 + |A(r)|(r3 cos^-t- 
Fssinc/i) -I- ihTiT2 — ifxT'iT^. The last two terms are 
the Zeeman field and the chemical potential, respec- 
tively. Choosing the representation for the matrices as 



Ti = az (E) ay, T2 = -cTz ® cTx, Ta = cr^ ® h, and 
Fs ~ —ay (g) I2 will then yield the original differential 
equations. The four matrices, together with Fq = CTz^CTz, 
are all Hermitian, and anticommuting among themselves, 
{Fi, Fj} = 25ij for z, J — 0, 1, 2, 3, 5. It is easy to see that 
in the absence of the last two terms the matrix Fq anti- 
commutes with the Hamiltonian, and thus ensures that 
the spectrum has the reflection symmetry around zero. 
In this limit there is a total of |ri| zero-modes for the 
vortex of vorticity n[ij, ^M,- It can be shown from the 
general properties of the Clifford algebra satisfied by the 
Dirac matrices that there also exists one antilinear op- 
erator that anticommutes with H even in the presence 
of the last two, ^ F1F2 and ^ F3F5 terms |35l|. In the 
representation defined above, in which Fi and F5 are the 
only two matrices among five which are imaginary, the 
antilinear operator in question is A = iTiT^K , where K 
denotes the operation of complex conjugation. Explicitly, 



A^ 



-lay 







K . 



(31) 



We recognize that the zero-mode ^0 hi Eq. p^ is its 
eigenvector with eigenvalue 1: A'^q = ^o- Here, it is 
worthwhile noting that Fq and A serve as the sublat- 
tice (chiral) and particle-hole symmetry operators, re- 
spectively. In the presence of finite chemical potential 
and/or Zeeman field, only the particle-hole symmetry is 
present and the JRD Hamiltonian belongs to the class D 
in the symmetry classifications. 36* In the case of vanish- 
ing chemical potential and zero magnetic field, the chiral 
symmetry is restored and the JRD Hamiltonian belongs 
to the class Dili. 

The effects of the Zeeman field on the zero- mode may 
be summarized as in Fig. HJa). For \h\ > y/fi'^ + A^, 
specified by the empty circle, the zero-mode is not nor- 
malizable, and hence disappears from the spectrum. The 
full circle, on the other hand, denotes the end of oscil- 
lations in the zero-mode in Eq. ([T5|) and ([TBI) . Lastly, 
the presence of a single Fermi surface is the key differ- 
ence in comparison with the same problem of the vor- 
tex in a noncentrosymmetric superconductor, in which 
the spin-orbit coupling is just the chiral kinetic energy, 
while the ordinary kinetic energy enter the Hamiltonian 
as i{k^ /2m — /i)F3r5. It is easy to see that while the 
operator A still anticommutes with the Hamiltonian, no 
zero-mode exists in the case for the s-wave symmetry of 
the order parameter 0. 

In the triplet case, the corresponding BdG Hamilto- 
nian has a momentum-dependent mass term. In the 
Dirac notation, the Hamiltonian is now H = Hjj + 
iMHd + ihTiT2 — i/iFsFs, where Ho and M are of the 
same form as, respectively, the kinetic energy and s-wave 
mass in the previous paragraph. It can be shown that 
the operator A still facilitates the reflection symmetry of 
the energy spectrum, due to the factor of i in the order- 
parameter term, which on the other hand is necessary to 
make it Hermitian. We showed here that the zero-mode 




(b) 
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FIG. 2: (Color online) The phases of zero-mode for h^ versus 
^^ in the (a) s-wave pairing and (b) p-wave pairing cases. 
The left-handed sides with respect to the empty circles denote 
the existence of the zero-mode, while the right-handed sides, 
including the empty circle, indicate that the zero-mode is not 
normahzable and hence does not exist. The blue regions stand 
for the fact that the wavefunction has, in addition to the 
decaying, an asymptotic oscillation part, which is gone when 
entering the green regions. 



does not exist when /i = 0, which can also be understood 
in terms of the symmetry of the zero-mode ^0 in Eq- 
(|24p under A: for positive and negative chemical poten- 
tial /I, we find that A^q = ±^0j respectively. In fact, 
the disappearance of the zero-mode at /z = suggests a 
topological transition[33], because the spectrum in Eq. 
^ in the translationally invariant case becomes gapless 
right at /J, = 0. In Fig. [2Jb) we summarize the effects 
of Zeeman field on the zero- mode. For h^ > fi^, there 
is no zero-mode, as again denoted by the empty circle. 
Similarly, the full circle at /i^ = /j,^/(l -t- v'^/vp), repre- 
sents the disappearance of oscillations under the overall 
exponential decaying in the zero-mode. 

In a real type-II superconductor, the Zeeman field h 
is not uniform and is proportional to the distribution 
of magnetic field B{r), a function riding on the scale 
of penetration length Ip. The obtained zero-mode ^0 
possesses two length scales kp^= v p / \/\iJ^^^-hF\ and 
^ — vp/A, which are assumed to be much smaller than 
the penetration length in the present problem. Therefore, 
it is valid to regard h(r) as a quasi-constant throughout 
the space. In fact, the distribution B{r) = ■^^KQ{r/X) 

for an isolated vortex, where (jjo is the magnetic flux. 
Hence, it could happen that near the core \h\ > |/i| as a 
result of \B\ ~ In— when r <C Zp[3J|. Correspondingly, 
the zero-mode will behave monotonically near the core, 
which is just the same as the the modified Bessel function 
with a length scale of VF/\h\. Nevertheless, the zero- 
mode will still be sustainable with the asymptotic h{r ^f 
00) within the range in Fig 121 

The mixed s-wave and p-wave symmetry is of relevance 
when it comes to the 2D surface considered here. In this 




FIG. 3: (Color online) Different phases for the zero-mode 
in the mixed pairing symmetry, as a function of fi and the 
parameter A_, assuming A+ is positive. The dotted axis is 
to exclude the special cases with fi — 0. The grey and blue 
regions denote different eigenvalues of A for the zero-mode 
^0 as specified in text. At the dashed line fi < 0, A- = 0, in 
between the two phases where A_ changes sign (as denoted 
by the green arrows) the localized zero-mode does not exist. 



case, the corresponding mass term is M + iMHe,, which 
still anticommutes with A and the principles outlined 
above can be applied. The amplitudes of the respective 
pairings are determined by their asymptotic values Aoo 
and UA- First, note that the Fermi surfaces for /x > 
and /i < have opposite helicities, and the corresponding 
zero-mode wave functions in Eq. (|29|) have the decaying 
length as vf/^+ and Wi?/|A_|, respectively. As shown 
in Fig. [31 the eigenvalue A of the zero-mode, A'i'o ~ 
A^o, may be different for opposite signs of the product 
A+A„. Assuming /j, is nonzero and A+ is positive, A 
remains the same when one changes /i = |/i| to /i = — |/i| 
if A_ > 0, which is similar to the s-wave case. However, 
A changes sign if A_ < 0, and this resembles the situation 
in the p-wave case. Besides, the zero-mode stops being 
normalizable when /i < and the parameters Aqo and 
kpVA have A_ vanishing, which corresponds to a gapless 
condition in Eq. (j31). This is specified by the dashed line 
in Fig. [21 The meaning of the parameters A± is clear 



when one considers the helicity basis, {a 



k4-'"k' 



}, which 



can be transformed back to the ordinary basis, {a 
with the following, 






}, 



''k+ 

4- 



V2 



le 



-i0k 




(32) 



where the angle 0k denotes the direction associated with 
the momentum on the plane. [llj, [32] Therefore, the time- 
reversal invariant pairing Hamiltonian can also be writ- 
ten as.llli 



k 



-i0k„t t 

"k+"-k-f 



-iA_e"'"^al_a\_) +h.c., 

(33) 

where the two independent order parameters A+._ are 
real, and corresponding to the pairing order parameters 
on the two helicity bands. Transforming back to the 
usual spin basis, the explicit singlet/triplet order param- 
eters are then given by, As(p) = (A+ ± A_)/2, respec- 
tively, which then yields the relation in Eq. p8|. It is 
worthwhile comparing with the noncentrosymmetric su- 
perconductors in which the presence of quadratic disper- 
sion term fc^/2m in the normal state leads to appear- 
ance of both helicity bands at the Fermi level. It follows 
that those with same (opposite) signs of A± resemble the 
singlet s-wave (triplet p- wave) superconductors. [11] Only 
the triplet superconductors are associated with topolog- 
ical states such as Andreev bound or zero-energy vortex 
bound states, which can exist near the topological de- 
fects. For superconducting surface Dirac fermions consid- 
ered here, the appearance of a single helicity band near 
the Fermi level is the key difference. [3^ All the pairing 
symmetries that we consider possess topological states 
as long as the superconducting gap is not closed by the 
Zeeman field or the decreasing chemical potential. Nev- 
ertheless, the mixed-parity pairing symmetry cases with 
same (opposite) signs of A± still resemble the pure s- 
wave (p-wave) case, manifested by the eigenvalue of an- 
tilinear spectrum-reflecting operator A as in FigO 

In conclusion, we determined the zero-modes for a vari- 
ety of the Jackiw-Rossi-Dirac Hamiltonians representing 
pairing of electrons on the surface of a topological in- 
sulator in presence of superconducting gaps of different 
symmetries. 
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